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Like a fairy-tale princess, trajectories around black holes can be sensitive to small disturbances, 
disturbance can lead to erratic orbits and an increased production of gravitational waves. 



We describe how a small 
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The high degree of symmetry found in the spacetime 
of an isolated black hole leads to regular geodesic mo- 
tion for orbiting bodies. Working on the assumption 
that small disturbances generally have small effects, it 
is often tacitly assumed that this idealised, textbook pic- 
ture carries over to real astrophysical situations. Here we 
want to emphasise that the idealised picture is not sta- 
ble against small perturbations, as the non-linearity of 
Einstein's equations tends to amplify small disturbances. 

The observation that small perturbations of an ide- 
alised black hole spacetime can lead to qualitative 
changes in the dynamics is not new. It has previously 
been noted that a range of perturbations lead to chaotic 
dynamics. The perturbations considered include addi- 
tional mass concentrations , magnetic fields || , 
gravitational waves [j?J and spin-orbit coupling || . Unlike 
the Kepler problem of Newtonian mechanics, essentially 
any perturbation of an isolated black hole spacetime will 
lead to chaotic orbits. This is because even the most pris- 
tine black hole spacetime harbours the seeds of chaos in 
the form of isolated unstable orbits. A small perturba- 
tion causes these unstable orbits to break out and infest 
large regions of phase space. Note that experience with 
Newtonian systems is very misleading. For example, the 
Kepler problem has more integrals of motion than are 
needed for integrability. Keplerian systems are thus im- 
pervious to small perturbations. In contrast, black hole 
spacetimes are at the edge of chaos, just waiting for the 
proverbial butterfly to flap its wings. 

Once it is realised that typical black hole - satellite sys- 
tems are chaotic in the strong field regime, we are forced 
to consider the consequences. One of the most immedi- 
ate consequences is that there will be no such thing as 
the "last stable orbit" ||. The boundary between sta- 
ble and unstable orbits will be fractal, and there may 
be large fractal tendrils of unstable orbits invading what 
would have been stable territory in an ideal black hole 
spacetime. This feature will be important in determin- 
ing when a binary system switches from inspiral to cata- 
clysmic collapse H . Another important consequence will 
be the increased production of gravitational waves due 
to the erratic motion of chaotic orbits. 

In what follows we will illustrate both of these effects 
in a simple model system. While our model does not de- 
scribe a real astrophysical situation, it does capture many 



of the salient features we expect to find in a relativistic 
binary system. 

For the purpose of illustration we will study the orbits 
of non-rotating satellites around an extreme Riessner- 
Nordstrom black hole. Almost identical results hold for 
motion around a Schwarzschild black hole perturbed by 
an orbiting third body Extreme black holes have 
the added advantage of allowing an exact generalisation 
to spacetimes with N extremal masses JIO|| . Rotating 
black holes bring with them a host of new instabilities 
which only amplify the points we wish to make. Simi- 
larly, rotating satellites introduce additional instabilities 
through spin-orbit || and spin-spin couplings. So while 
our model is chosen on the grounds of simplicity, it is 
likely that more realistic models will be even more sensi- 
tive to small perturbations. Moreover, a central feature 
of non-linear dynamics is universality: the details of the 
dynamics are less important than the general structure 
of the phase space trajectories. For example, a stochastic 
layer will lead to an increase in the gravitational wave lu- 
minosity regardless of what lead to the formation of the 
stochastic layer. Because of this, the effects we describe 
for our particular model will enjoy wider applicability. 

Our unperturbed spacetime is described by the metric 
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and electromagnetic potential 
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is the metric of a 2-sphere. 



Here dVt 2 = d9 2 + sin 2 
Into this spacetime we introduce a spinless satellite of 
mass [i and consider its orbits. To keep our model simple 
we will neglect the spacetime curvature caused by the 
satellite. The rotational symmetry allows us to restrict 
our attention to orbits in the plane 9 — tt/2. Moreover, 
invariance under time translations and spatial rotations 
leads to conservation of the satellite's energy [iE and 
angular momentum fiL. This allows the motion to be 
reduced to the radial equation 
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where r is the proper time along the trajectory and 
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We will be concerning ourselves with bound orbits {E < 
1), but similar conclusions hold for unbound orbits. In 
Fig. 1 the effective potential V(r) is displayed for three 
values of L. The energy E = 0.94 is also displayed. We 
see that some orbits can reach the event horizon at r = 
while others are protected by an angular momentum bar- 
rier. For E = 0.94 this occurs at the critical angular 
momentum L c — 2.982Af. A trajectory with these pa- 
rameter values can follow an unstable circular orbit at 
r c = (L 2 - 2M 2 - L cy /L 2 c - 8M 2 )/2M. As emphasised 
by many authors, it is the existence of such unstable pe- 
riodic orbits that make black hole spacetimes sensitive to 
small perturbations. 
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leading to 
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The new metric is an exact solution to the Einstein- 
Maxwell field equations belonging to the Majumdar- 
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10H family of spacetimes. The leading cor- 
le extreme Riessner-Nordstrom metric are 



AU = - l + -cos< 



+ ^2 [3 cos 2 0-1] 
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The dipole term, cos </>, can be removed by a coordinate 
transformation, while the quadrupole term, (3 cos 2 4> ~ 
1), is analogous to a perturbation of the Schwarzschild 
metric studied by Moeckel || . He described how a small 
body orbiting a Schwarzschild black hole causes the inner 
orbits to become chaotic. 
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FIG. 1. The effective potential V(r) for different values of 
the angular momentum. In order of decreasing height, the po- 
tentials correspond to L = 3.1M, 3M, and 2.9M respectively. 
The solid line is the fixed energy E = 0.94. 
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FIG. 2. Basin boundaries in the (L, r) plane. 



Previous studies have used local methods to assess the 
onset of instability. Here we apply a more revealing 
global approach based on stability basins. Such global 
techniques are increasingly being used in engineering ap- 
plications in order to survey the full range of stable and 
unstable configurations of circuits, bridges and boats |pl| . 
When systems are chaotic, stable regions of phase space 
can be invaded by chaotic tendrils. These tendrils exhibit 
a complicated fractal structure. 

To illustrate how small perturbations can lead to im- 
portant changes in the dynamics we introduce a small 
"pea" with mass m and charge q — m at a distance r = R 
from the black hole. The new metric is obtained by re- 
placing 1 + M/r by 



Since our perturbed metric is static, the satellite's en- 
ergy remains a conserved quantity. In addition, the mod- 
ified potential is independent of 6, so pe is conserved. 
This means we can continue to study trajectories in the 
plane 9 — it/ 2 without loss of generality. However, the 
4> dependence of the potential breaks rotational invari- 
ance so the satellite's angular momentum is no longer 
conserved. A trajectory with L(0) > L c can now evolve 
to one with L(r) < L c , thus allowing the possibility of 
capture by the black hole. Conversely, a trajectory that 
was destined for captured might now gain enough angu- 
lar momentum to avoid capture. For chaotic trajectories 
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this gain or loss of angular momentum can depend sensi- 
tively on initial conditions, leading to complicated fractal 
boundaries separating the stable and unstable outcomes. 

Here we are working under the usual assumption that 
we can use the dynamics of the dissipationless system 
to predict which orbits will be unstable when gravita- 
tional radiation is included as a form of dissipation. This 
amounts to a kind of adiabatic approximation which is 
valid when the energy loss per orbit is much less than the 
energy of the satellite. 

In Fig. 2 we display the stability basin in the (L, r) 
plane for a particle with energy E = 0.94. A small pea 
with mass m = M/100 has been introduced at r = R = 
10M. Particles are started from an initial position r, with 
"angular momentum" L — r 2 U 2 d<f)/dT. The trajectory 
is then evolved numerically. If the particle is captured 
by the black hole we colour the initial position white. 
If the particle achieves a stable orbit (defined here to 
be 100+ orbits without being captured), then we colour 
the initial position black. There are two major features 
of note in Fig. 2. The first is the irregular boundaries 
between stable and unstable outcomes. The second is the 
broad swath of unstable trajectories between L = 3M 
and L = 3.15Af. Before the pea was introduced, all 
trajectories with L > L c = 2.982M were stable. 
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FIG. 3. A detail of Fig. 2 near L = L c . 

The stability basins for the unperturbed spacetime are 
separated by a smooth line at L — L c . In Fig. 3 we 
display a detail of Fig. 2 which clearly shows that the 
smooth boundary near L c has been replaced by a com- 
plicated fractal structure. Note that most orbits with 
L < 2.995M have been rendered unstable and thrown 
into the black hole. In Fig. 4 we see that the central 
basin of instability also has a fractal boundary. Fractal 
structures provide a gaug e invariant signal of chaos in 
general relativity fl ^2] , ^3[ . 



But aside from producing nice pictures, what does 
chaos add to the physics of satellites orbiting black holes? 
One interesting possibility is an enhancement of the grav- 
itational wave output. Chaotic trajectories will tend to 
produce gravitational waves with increased luminosity 
and amplitude, and erratic variations in amplitude. The 
power radiated will exceed that of stable orbits since, in 
loose terms, the luminosity is proportional to the rate 
of change of acceleration. Chaotic orbits often exhibit 
highly variable accelerations while nearly circular orbits 
have gently varying acceleration profiles. In fact, even 
elliptical orbits produce considerably more gravitational 
radiation than circular orbits due to the sharp turn at 
their point of closest approach. Unfortunately, most or- 
bits are expected to be nearly circular by the time t hey 
reach their swan song in the LIGO detection band p4[ . 
This is because tidal friction and gravitational radiation 
act to circularise the orbit. To get more powerful grav- 
itational wave signals we need something to destabilise 
this picture. Enter chaotic resonances. 
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FIG. 4. A detail of Fig. 2 showing the boundary of the 
band of instability. 

As a realistic binary system spirals inwards, various 
relativistic instabilities start to become important. The 
instabilities might be caused by spin-orbit or spin-spin 
coupling, or by external mass distributions such as the 
third body studied here. While these resonances are typ- 
ically restricted to isolated bands in phase space, the in- 
spiralling satellite is likely to run across an unstable band 
as its energy and angular momentum are reduced by the 
emission of gravitational waves. The resulting chaotic or- 
bit would then provide a boost to the gravitational wave 
output. 

To get a feel for this chaotic enhancement we can 
calculate the wave amplitude and power radiated using 
the quadrupole approximation. Since the orbits we are 
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studying are fairly relativistic, the quadrupole approxi- 
mation can only provide a qualitative picture. The power 
radiated by the satellite is given by |Q 

P=\<{4h 3 f>, (8) 

where <> denotes the average over several orbits and 7^ 
is the reduced quadrupole moment of the satellite's orbit 
(i,j = 1,2,3): 

Jy = (xiX 3 - ^r 2 )fi . (9) 

The notation d 3 is shorthand for d 3 /dt 3 . The direction- 
averaged wave amplitude is given by 

A oc ^(dtlijMlv) . (10) 

Switching to cartesian coordinates, and considering or- 
bits in the plane z = (9 = ir/2) we find 

A cc (HdUxy)f + [d 2 t x 2 - d\y 2 f + {d 2 t x 2 ){d 2 t y 2 )) 1/2 . 

In Fig. 5 we display the amplitude of three represen- 
tative trajectories with E = 0.94, normalised against 
the stable circular orbit of the unperturbed system with 
L = 3.1518M. 
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FIG. 5. The amplitude of circular (dashed line), elliptic 
(dotted curve) and mildly chaotic (solid curve) satellite orbits. 

The dotted line shows the variation in amplitude of a 
precessing elliptical orbit of the unperturbed system with 
L = 3.117M. The solid line is for an orbit with the 
same initial conditions as the regular elliptic orbit, but 
this time for the perturbed system. The variation in 
amplitude is increased by up to a factor of two over the 
unperturbed system. The elliptic and chaotic orbits have 
power outputs 1.48 and 2.01 times larger than the stable 
circular orbit. This implies that the small perturbation 
caused by the "pea" increased the power output by 36%. 



To better assess the significance of this result, we need 
to understand the types of trajectories that give rise to 
increased power output, how generic these might be in 
physically reasonable spacetimes and how large the ef- 
fect can be. To this end, we display in Fig. 6 a Poincare 
section for trajectories with E — 0.94. The Poincare sec- 
tion reveals a combination of unbroken KAM tori, island 
chains, cantori and thin stochastic layers. These features 
are typical for mildly chaotic systems. The chaotic orbit 
shown in Fig. 5 forms the thin stochastic layer that makes 
up the outermost ring of the Poincare section. While this 
is orbit is amongst the most chaotic we found, it only er- 
godically wanders over a very small band in phase space. 
The unbroken KAM tori prevent orbits from becoming 
highly erratic. This in turn limits the gravitational wave 
output. Moreover, orbits of the perturbed system that lie 
on the unbroken KAM tori experience much smaller in- 
creases in power output. For example, the outermost un- 
broken KAM tori in Fig. 6 produces just 5% more power 
than the analogous orbit of the unperturbed spacetime. 
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FIG. 6. A Poincare section for trajectories with E — 0.94. 

These observations suggest that the power output is 
greatest for orbits that fill the largest regions of phase 
space. In the model spacetime we have been studying 
the dynamics is dominated by unbroken KAM tori, and 
this limits the size of any stochastic layers. In addition, 
we studied a 4-dimensional Hamiltonian system so the 
unbroken KAM tori partition phase space. Systems with 
fewer KAM tori, higher dimensional dynamics, and/or 
dissipation typically have trajectories that wander over 
large regions of phase space. In such systems the ef- 
fects we have been describing will be larger and more 
widespread. 

The literature already contains several realistic black 
hole models where highly erratic orbits have been found. 
For example, some of the trajectories shown in Fig. 1 of 
Ref. ||, and Fig. 4(f) of Ref. || form wide stochastic 
layers. We are currently modelling these spacetimes to 
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see how large the gravitational wave enhancement can 
be. 

Using a simple model we have shown that chaotic insta- 
bilities will affect a satellite's transition from inspiral to 
plunge, and cause an increase in gravitational wave pro- 
duction. Both of these features could be important when 
producing gravitational wave templates for the LIGO and 
LISA detectors. 

We thank Sam Drake, Eric Poisson and Janna Levin 
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the Australian Research Council. 
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